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We show that the recently demonstrated technique for generating stationary pulses of light [Nature 
426, 638 (2003)] can be extended to localize optical pulses in all three spatial dimensions in a 
resonant atomic medium. This method can be used to dramatically enhance the nonlinear interaction 
between weak optical pulses. In particular, we show that an efficient Kerr-like interaction between 
two pulses can be implemented as a sequence of several purely linear optical processes. The resulting 
process may enable coherent interactions between single photon pulses. 



Techniques that could facilitate controlled nonlinear 
interactions between few-photon light pulses are now ac- 
tively explored 0]. Although research into fundamental 
limits of nonlinear optics has been carried out over the 
last three decades, there is renewed interest in these prob- 
lems in part due to e.g. potential applications in quantum 
information science |2( . In general, such interactions be- 
tween few-photon pulses are difficult to achieve, as they 
require a combination of large nonlinearity, low photon 
loss and tight confinement of the light beams U . In ad- 
dition, long atom-photon interaction times are required. 
Simultaneous implementation of all of these requirements 
is by now only feasible in the context of cavity QED 0. 

In this Letter we describe a novel method for achiev- 
ing nonlinear interaction between weak light pulses. Our 
method is based on a recently demonstrated technique 
in which light propagating in a medium of Rb atoms 
was converted into an excitation with localized, station- 
ary electromagnetic energy, which could be held and re- 
leased after a controllable interval. This is achieved by 
using Electromagnetically Induced Transparency (EIT) 
to coherently control the pulse propagation. We show 
here that this method can be extended to confine station- 
ary pulses in all three spatial dimensions. This, in turn, 
can be used to strongly enhance the nonlinear interac- 
tion between weak pulses of light. Specifically we demon- 
strate that an efficient Kerr-like interaction between two 
pulses can be implemented as a sequence of linear opti- 
cal processes and atomic state manipulations. Coherent, 
controlled nonlinear processes at optical energies corre- 
sponding to a single light quanta appear feasible. 

Before proceeding, we note that the present work 
is closely related to recent studies on the resonant 
enhancement of nonlinear optical phenomena via EIT 
0, 0, 0, O]. The essence of these studies is to uti- 
lize steep atomic dispersion associated with narrow EIT 
resonances. In such a system, a small AC Stark shift 
associated with a weak off-resonant pulse of signal light, 
produces a large change in refractive index for a reso- 
nant probe pulse. In order to fully take advantage of 
this process, long interaction times between signal and 
probe pulses must be ensured. Although the latter can 
be achieved by reducing the group velocities of two in- 
teracting pulses by equal amounts in practice this 



results only in a modest increase of the nonlinear optical 
efficiency since reduction of the group velocity is accom- 
panied by a corresponding decrease of the light energy in 
the propagating pulse. Moreover, such nonlinear inter- 
action is accompanied by pulse distortion, which poses a 
fundamental limit to nonlinear interactions. In contrast, 
the technique presented here allows for long interaction 
times (associated with stationary light pulses) without 
proportional reduction of the photon energy. Further- 
more, our technique enables light localization in all three 
spatial dimensions in the presence of atoms, and allows to 
entirely avoid competing effects such as pulse distortion. 

The ideas discussed in Refs. allow one to localize 
and hold a pulse of light within a stationary envelope 
along the propagation direction. In practice, focused 
pulses will undergo diffraction in the transverse direc- 
tions. In order to prevent diffraction, it is necessary to 
confine the signal beam in transverse directions. This can 
be achieved, e.g. by using a hollow core photonic crystal 
fiber filled with an active medium of resonant atoms |l3| . 
Such an approach is particularly attractive in that it en- 
sures a single-mode beam quality for interacting beams. 
The unwanted interactions of atoms with fiber walls can 
be avoided by using atom guiding techniques 0. Al- 
ternatively, the signal pulse guiding can be accomplished 
by using focused control beams. Shaped control beams 
can be used to create a transverse variation of the index 
of refraction, enabling waveguiding and confinement of 
light pulses to small transverse dimensions. 

To be specific we consider a medium of length L con- 
sisting of an ensemble of N three-level atoms in the A 
configuration, with two metastable lower states, as shown 
in Fig. la. The ground states |<?),|s) are coupled to 
the excited state |e) via a control field applied on res- 
onance with the \s) — > |e) transition and a weak quan- 
tized signal field close to resonance with the \g) — > |e) 
transition. The control field consists of two counter- 
propagating fields with spatially and temporally vary- 
ing Rabi frequencies Q± (r, t) , so that the control field 
Rabi frequency is Q(r,t) = VL + (r,t)e lkcZ + VL-(r,t)e~ lkcZ , 
where k c = n c u> es /c, n c being the background index of 
refraction at the frequency co es . The corresponding sig- 
nal fields have slowly varying envelopes E±(r,t), so that 
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FIG. 1: a. Three-level atoms in A-configuration, with auxil- 
iary levels \s'), |e'). Forward and backward propagating con- 
trol fields with Rabi frequencies fi±, and weak signal field £±. 
b. Three-dimensional confinement and waveguiding of signal 
light due to transverse intensity profile and longitudinal mod- 
ulation of control field. 



the signal field is 
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where k s = n s Lu eg /c, V is the quantization volume, and 
n s is the background refractive index at the frequency u> eg 
due to off-resonant atomic levels. Note that in practice 
n s can be tuned, e.g. by changing the light polarization. 

We describe the atomic pro per ties with slowly 
varying collective operators [15j ct m „ (r , t) = 

W?Ef=i\^)jW\^ iU, " vt wh ere uv = (E^ - E v )/h, 
and where Np is the number of atoms in a small but 
macroscopic volume around position r. We define the po- 
larization operator to be P(f,t) = \fNa ge (r,t), and the 
spin flip operator S (r , t) — y/~N a gs (r, t) . In the present 
situation of weak signal fields and strong control fields, 
most atoms are in state \g), with a few spin-flipped atoms 
in |s), so that the polarization and spin flip operators 
obey bosonic commutation relations 0|. Associated 
with the forward/backward propagating fields are slowly 
varying polarization envelopes P± (r , i) , so that the total 
polarization is P(f, t) = P + (f, t)e lksZ + P-(f f , t)e~ tksZ . 

Letting the wavevector mismatch between co- 
propagating signal and control fields be AK = k s — k c = 
(n s ui eg — n c uj es )/c, and defining £± = E±e ±lAKz and 
V± = P±e ±lAKz , the Heisenberg equations of motion for 
the slowly varying operators £± (f, t) can be written 
as (in the paraxial approximation) 
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£± = iAKc£± + igVNV±,(2) 
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■j^2 is the transverse Laplacian, and 
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is the atom-field coupling con- 



where Vy = V 2 

where g — p ( 

stant, p being the dipole matrix element for the \g) — |e) 
transition, and V the quantization volume. 

Following 0,^3 we introduce two components 4>± of 
a coupled excitation of light and an atomic spin wave 



("dark-state polariton") corresponding to forward and 
backward signal fields respectively. In the experimentally 
relevant case of small group velocities 0, 0] the polari- 
ton components are represented by ^± = g^/N£±/Q±. 
In the adiabatic limit of slowly varying pulses, disregard- 
ing the slow decay of ground state coherence, and Fourier 
transforming (dt — > —ioj), we find 
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and £ = 3-£L - We have also assumed that k Q a 3> 1, 
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where a is the typical transverse size of the control beams. 
These equations describe two slow waves that are cou- 
pled due to periodic modulation of atomic absorption 
and group velocity. The terms containing £ on the right 
hand side of Eqns. © are proportional to the absorption 
coefficient £ near resonant line center. F±(r,ui) are noise 
forces associated with dissipation. When £ is large these 
terms give rise to the pulse matching phenomenon 0: 
whenever one of the fields is created the other will adjust 
itself within a short propagation distance to match its 
amplitude such that * + - #_ — > @]. 

In order to achieve transverse confinement of light 
pulses, we take into account the transverse dependence 
of the control field intensity and the resulting variation of 
the index of refraction. For a focused control beam, the 
intensity decreases with distance from the optical axis, so 
that for negative two-photon detuning (as is necessary for 
phasematching), the index of refraction decreases with 
distance from the optical axis. This leads to waveguid- 
ing of the signal light. The combination of waveguiding 
with strong coupling of forward and backward propagat- 
ing modes, permits the complete three-dimensional con- 
finement of light pulses in the medium. 

We assume a transverse spatial variation of the control 
field intensity, e.g. for a weakly focused gaussian beam, 
\n±(r)\ 2 = |O ± (0)| 2 e~^/ Q ) 2 . Expanding for r « a, we 
have T](r) = rjo[l + (r/a) 2 + •••]. We consider trial 
solutions of © of the from = A+e^"'^ I and 
= A_ e -^-(r/-R) 2 where /3, R and A + /A_ are deter- 
mined by requiring that the coefficients of different pow- 
ers of r vanish independently. With these requirements 
we find the two solutions R = oo and 
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The eigenvector of the finite solution has A+/A^ = 1, 
which corresponds to the stationary pulses |6j for which 
ty + — — * 0. We also find the dispersion relation 
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AKc^j - i^- + (a+ - a-)cp. (5) 



In the time-domain this corresponds to propagation at a 
group velocity v g — c^—- 



that can be controlled by 
adjusting the intensities I±(i) oc |ri±(i)| 2 of the counter- 
propagating control fields. Due to the imaginary term, 
there is also a slow spreading of the stationary pulse at a 
rate 81/1 ~ \J c 2 t / (rj^l 2 ) , which determines the maximal 
trapping time of the stationary excitation. 

We are interested in a simultaneous solution of Eqs.© 
when a+ ps a_ and the optical depth £ is large. 
This yields the radius of the stationary pulse R = 

1/2 

02 1 / 4 + AKkoa 2 — l] , which under conditions 
of strong confinements AKkoa 2 ^S> 1 results in 



R = a[2/(AKk a 2 )} 1/4 . 
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Hence, in an optically dense medium (£L/c 3> 1) a sta- 
tionary excitation confined in all three dimensions can 
be controllably created. To be specific, for atomic Rb 
(A = 0.8^m) at density n — 10 14 cm~ 3 , we take the back- 
ground refractive index due to off-resonant levels to be 
n — 1 = 1.2 • 10~ 2 18]. For a Gaussian control beam 
with waist a ~ 100/im at the center of the atomic cell 
( the corresponding Rayleigh range is zq = 3.9cm), we 
find that the guided mode radius is R — 13/im (for 
which the Rayleigh range is zq — 0.06cm), so that the 
diffraction-free range is extended by a factor of 60. When 
the control beam is chosen in the form of a non-diffracting 
Bessel beam [I9j, waveguiding over much longer propa- 
gation distances is possible. For example, with a beam 
in which the radius of the first lobe of the Bessel function 
is a = 20/im, the guided mode radius is R ~ 5.7/im. In 
practice, this allows confinement over tens of cm, whereas 
in free space the corresponding Rayleigh range would 
only be 0.01cm. Reducing the control beam radius until 
R ~ a, gives an estimate of the smallest guided mode 
achievable, which for an index of n — 1 = 1.2 • 10~ 2 is 
Rmin = 1-6/xm. 

We next turn to the nonlinear optical interaction be- 
tween two weak pulses of light. A notable feature of 
the step by step process described below, is that it con- 
sists of a sequence of purely linear optical interactions 
and atomic state manipulations, leading to an effective 
optical nonlinear interaction. As shown in the timing 
diagram Fig. 2a, a light pulse (signal pulse) travelling 
through the atomic medium is initially stored in the \g)(s\ 
coherence. Next, a Raman or microwave n pulse (RA) 
transfers the population from \s) to \s') (see Fig. 2a), 
thereby transferring the stored excitation to the spin ex- 
citation S' — ^/N&gs'. The signal pulse is stored as a 



spin wave a gs i (r, t) = U s (p)S' (z,t)/VN, where U s (p) de- 
scribes the transverse spin wave mode. 

A second light pulse (probe pulse) is then sent through 
the medium and stored in the coherence S. The spin exci- 
tation associated with the probe pulse is then converted 
into a stationary excitation. The latter is then moved 
through the stored spin excitation associated with the 
signal pulse. Probe light in the modes £± interacts dis- 
persively with atoms in level |s') (see Fig. la), thereby 
acquiring a phase shift proportional to the number of 

excitations S'^S', leading to an effective Kerr-type non- 
linearity 

We focus on the situation where the transverse spa- 
tial size of the signal spin wave S' is much smaller 
than the probe transverse size. Under these assumptions 
and writing the guided mode transverse dependence as 



y±(r,z,t) = e 
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ip±{ z -,~t)i we find, in terms of the 



polariton components tp± 1 

(d t + cd z )ijj + = -r)d t (a + 4i + + «_?/>_) - a_£(V> + - ?/>-) 
+ if3S' Jf S'iJ> + + F + (z,t) (7a) 
(d t - cd z )ij>- = -r)d t (a + 4! + + «_?/>_) + a + £(4> + - 

(7b) 
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where f3 = ^-(1 + »7 / 'A), with g = g^jp, A is the quan- 
tization area, and A is the detuning of the fields £± from 
the optical transition \s') — > |e') (see Fig. 2a). The in- 
teraction of the localized excitation S' with the guided 
modes ip± does not depend on the transverse coordinate, 
and the effective transverse area corresponds to the mode 
area 7ri? 2 , in complete analogy to the interaction of lo- 
calized atoms with the field mode in cavity QED jij. 

In the stationary pulse configuration [(|, the station- 
ary excitation is bound to the spin wave and S ~ 
— (^aj r '4> + +a_^_^. Solving adiabaticaly, in the limit 

of large optical depth £ and writing v g = (a + — a-)c/r/, 
we find 



[d t + v g d z ]S = i-[S ,i S']S 



V 

d t S' = i-[&S\S' 



(cd z 



S (8a) 
(8b) 



where we have ignored for now absorption and the asso- 
ciated noise. 

To solve for S(z,t) and S'(z,t), we first ignore pulse 

spreading. Let hx(z) = S'^S' (independent of t), and 
fi2(z, t) — S^S (which depends only on the variable t' — 



4 



a. 

E 
S 
S' 
ED 



b. 



BD 
RA 



"L 



_J A^. 



FIG. 2: Nonlinear optical interaction with weak pulses as 
sequence of linear operations, a. Timing diagram: forward 
(FD), backward (BD), and Raman or microwave (RA) inten- 
sities vs. time. b. Illustration of Kerr interaction between 
slowly propagating stationary pulse S and stored excitation 
S' leading to phase shift (represented as changing color) of 
pulses. 



z/vg). We find 



S{z, t) — exp 



S'(z, t) = exp 
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V v g J 
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S(z-v g t,0) (9a) 



S'{z,0) (9b) 



When the slowly moving pulse S has completely tra- 
versed the stored spin coherence S', the phase shift is 



<t>s 



Rc[/3] 



LNgi (where N$i is the number of excitations 



initially stored in S'). The phase shift is proportional 
to the interaction time, i.e. inversely proportional to the 
group velocity of the slowly moving pulse S. To estimate 
the maximal phase shift, we note that the group veloc- 
ity must be large enough that v g t > l s i, where l s i is the 
length of spin coherence envelope. Also, non-adiabatic 
corrections due to the pulse spreading term in QSaj l should 

be small, so that \.c t ^ 1. Putting these two condi- 
tions together yields 



4>s < di 
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where the resonant scattering cross-section is a — -p-A 2 . 
Note that the nonlinear phase shift scales linearly with 
the optical depth do, in contrast to scaling with V^o for 
the case of two slowly propagating pulses , in which 
case pulse distortion effects are also significant. 

Specifically, a 300^m long, cigar shaped cloud of cold 
87 Rb atoms confined in an optical dipole trap at a density 
of n ~ 10 14 cm -3 has an optical depth in excess of do ~ 
10 3 . Similar optical depth can be potentially achieved by 
guiding cold atom clouds of smaller density in a photonic 
crystal fiber [lil IHj . Taking the guided mode radius 
R ~ 2/im and accounting for absorption losses, we choose 



the detuning to be A ~ I67 and find that a phase shift 
of 05 ~ 7r is achievable due to a single stored quantum. 
Under these conditions the two-photon loss probability 
is a few percent. 

To summarize, we have shown that three-dimensional 
confinement of light pulses is possible by combining the 
technique of stationary light pulses with the transverse 
light guiding. This technique can be used to engineer 
efficient nonlinear optical interactions leading to signif- 
icant phase shifts for weak optical pulses. Such inter- 
actions have interesting applications ranging from QND 
measurements j'2(l of few-photon pulses to quantum in- 
formation processing 0. 
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